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ABSTRACT 



Aims. We aim at investigating the formation of jet-like features in the lower solar atmosphere, e.g. chromosphere and transition region, 
as a result of magnetic reconnection. 

Methods. Magnetic reconnection as occurring at chromospheric and transition regions densities and triggered by magnetic flux emer- 
gence is studied using a 2.5D MHD code. The initial atmosphere is static and isothermal, with a temperature of 2 x 10"* K. The 
initial magnetic field is uniform and vertical. Two physical environments with different magnetic field strength (25 G and 50 G) are 
presented. In each case, two sub-cases are discussed, where the environments have different initial mass density. 
Results. In the case where we have a weaker magnetic field (25 G) and higher plasma density (N^ = 2 x 10" cm"''), valid for the typ- 
ical quiet Sun chromosphere, a plasma jet would be observed with a temperature of 2-3 x 10^ K and a velocity as high as 40 km s"'. 
The opposite case of a medium with a lower electron density (A'^ = 2 x lO"' cm"'), i.e. more typical for the transition region, and a 
stronger magnetic field of 50 G, up-flows with line-of-sight velocities as high as ~90 km s"' and temperatures of 6 x 10^ K, i.e. upper 
transition region - low coronal temperatures, are produced. Only in the latter case, the low corona Fe ix 171 A shows a response in 
the jet which is comparable to the O v increase. 

Conclusions. The results show that magnetic reconnection can be an efficient mechanism to drive plasma outflows in the chromo- 
sphere and transition region. The model can reproduce characteristics, such as temperature and velocity for a range of jet features like 
a fibril, a spicule, an hot X-ray jet or a transition region jet by changing either the magnetic field strength or the electron density, i.e. 
where in the atmosphere the reconnection occurs. 

Key words. MHD - Sun: chromosphere - Sun: transition region - Sun: corona - Sun: magnetic fields 
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1. Introduction 

A variety of jet-like features such as spicules, fibrils, surges, 
Ellerman bombs, EUV/X-ray jets etc., as seen in various atmo- 
spheric regions, are observed in the solar atmosphere. Spicules 
are relatively thin, elongated jet-like structures best viewed at the 
solar limb in Ha or Ca II images as bright features against a dark 
background. The close correlation between observed properties 
of limb spicules and other on-disk features such as mottles, fib- 
rils and straws has prompted many authors to suggest that these 
may be counterparts of each other (Tsiropoula & Schmieder 
1997). As viewed in the optical, the classical spicule is observed 
to reach heights of 6500-15 000 km (Beckers 1968; Withbroe 
1983) with an average lifetime of 5 minutes and average plasma 
velocities of 25 km s"' . Recently, using high-resolution observa- 
tions in Ca II H (3968 A) from the Solar Optical Telescope (SOT) 
on Hinode, De Pontieu et al. (2007) suggested at least two types 
of spicules, the classical spicule and a more dynamic one. 

Tavabi et al. (2011) suggested four types of spicules based 
on their diameter, ranging from 0.3" (220 km), 0.5" (360 km), 
0.75" (550 km) to 1.15" (850 km). Typically, they show a suc- 
cession of upward and downward motions. The more dynamic 
spicule develop and disappear on timescales of 10-60 s, with 
velocities sometimes exceeding 100 km s"'. 

Spicules are also seen in UV and EUV lines (Dere et al. 
1989) and thus reach at least to transition region tempera- 
tures. Cooketal. (1984) using HRTS observations taken in 



C IV 1550 A showed that these structures show tilted features 
which was interpreted as rotational velocities of approximately 
50 km s"'. In a more recent paper, Madjarska et al. (2011) 
showed that this may be better interpreted as a multi-strand 
structure with up-flows and down-flows. 

Several authors have looked at magnetic reconnec- 
tion as a driving mechanism for hot jets (Jin et al. 1996; 
Innes&Tdth 1999; Galsgaard et al. 2005; Nishizuka et al. 
2008; Patsourakos et al. 2008; Murray et al. 2009; 
Rosdahl & Galsgaard 2010; Pariatetal. 2010, and refer- 
ences therein). Roussev et al. (2001a,b,c) performed 2D 
Magneto-hydrodynamic (MHD) simulations of transition region 
jets deriving blue-shifts of the order of 100 km s"'. Several 
initial physical environments were studied. However, the plasma 
/3 on the current sheet was the same in all the cases with the 
maximum velocity of the blue-shifted jets eventually reaching 
almost the same value in all cases, although they were different 
at the beginning of the experiments. 

Various authors have suggested that spicules can be 
driven by waves (magnetoacoustic or Alfven) (Hansteen et al. 
2006; De Pontieu et al. 2007; Heggland et al. 2007), while 
Sterling et al. (1993), Karpen et al. (1995), and Heggland et al. 
(2009) have considered magnetic reconnection as a plausible 
candidate for driving chromospheric jets. In all the models men- 
tioned above, spicules are produced but their velocities are small 
(~ 25 km s"'). Martfnez-Sykora et al. (2010) explored 3D simu- 
lations of flux emergence, and reported on more dynamic fea- 
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tures that reached coronal temperatures. In their model, large 
tangential discontinuities of magnetic field is necessary to cre- 
ate the conditions for accelerating the jets. These large magnetic 
field gradients lead to strong Lorentz forces, causing horizontal 
flows which reach a "wall" of strong vertical magnetic field. The 
wall squeezes the plasma and forces a large increase in the pres- 
sure, deflecting the horizontal flows and pushing them to move 
vertically. 

Table 1. Parameters of the four initial states, namely Al, A2, Bl 
and B2, including the characteristic values of the Alfven velocity 
(va) and plasma fi in each state. 



PbQ 




2 




1 






Al 




Bl 


1 


va = 


244 km s"' 




= 122 km 




P = 


1.1 X 10-2 


P 


= 4.4 X 10-' 






A2 




B2 


0.1 


Va = 


772 km 


Va 


= 386 km S-' 




p = 


1.1 X 10-3 


P 


= 4.4 X 10-3 



Pariat et al. (2009) proposed a 3D model for solar polar jets, 
where the magnetic twist was taken as the jet driver. The release 
of magnetic twist onto open field lines by magnetic reconnection 
resulted in high-speed jets. Their work reproduced helical struc- 
tures observed in some polar-hole jets (Patsourakos et al. 2008). 
Jiang et al. (201 1) examined the influence of different initial re- 
connection angles by using a 3D MHD model where thermal 
conduction and gravity were neglected. Fan-shaped jets moving 
along the magnetic guide field were obtained. 

Recently, Ding et al. (2010) used a 2.5-dimensional resis- 
tive MHD model in Cartesian coordinates to investigate mag- 
netic reconnection in the low atmosphere, e.g. chromosphere, 
discussing the implications for jet features at transition-region 



temperatures. They showed that faster and hotter outflows could 
be obtained if a physical environment with lower mass density 
was considered, however, the highest temperature of the plasma 
heated by magnetic reconnection is only 5x10^ K in their study. 

Here, we expand the work of Ding et al. (2010) varying the 
electron density and field strength as applied to a larger range of 
jet features, including fibrils, spicules, chromospheric jets and 
transition region jets. This work uses a larger grid model (see 
later) than the previous study. Here, we also look at jet for- 
mations as viewed in different parts of the lower solar atmo- 
sphere. We derive the response of three spectral lines with for- 
mation temperatures from the lower transition region to the low 
corona, as obtained in the foot-points and in the jet itself. The 
2.5-D resistive MHD model is briefly described in Section 2. 
Section 3 gives the numerical results and the derived line pro- 
files. Conclusions and discussion are drawn in Section 4. 

2. Physical Model and Numerical Methods 

2.1. Basic Equations 

A 2.5-dimensional resistive MHD model in Cartesian coordi- 
nates is used here. The MHD equations are the same as pre- 
sented in Ding et al. (2010), where p, v, tfr, B~, T, p, j, B, g, tj, j, 
Pd, Q, Lr are the mass density, flow velocity, magnetic flux func- 
tion, j-component of magnetic field, temperature, gas pressure, 
electric current density, magnetic field, gravitational accelera- 
tion, dimensionless magnetic diffusivity, adiabatic index, char- 
acteristic ratio of the gas pressure to the magnetic pressure de- 
duced from the basic units, heat conduction, and radiative losses, 
respectively. And, j, B, Q, Lr are explicitly expressed by 

p =pr, j = V xB, andB = V X (i/rz) + (1) 

g = V ■ [T^l\^ ■ vr)B/B-], and = p^A(r), (2) 

where A(r) is the radiative loss function. Here, optical depth 
effects are considered: At T < 2 x 10"* K, K{T) is reduced 
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(a) Case Al: 50 G, 2x10" cm"' 




(b) Case A2: 50 G, 2xl0'° cm"' 
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(c) Case Bl: 25 G, 2x10" cm"' 



(d) Case B2: 25 G, 2xl0'° cm" 
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Fig. 2. Evolution of the magnetic field {solid lines), temperature (colour), and velocity (arrows) at four times for case Al (a), A2 (b), 
Bl (c), and B2 (d). The characteristic parameters of magnetic field strength and electron number density are given above each case 
plot. Note that the colour-coding is different for each run. 



to zero; At 2 x 10'' < T < 10^ K, A(T) suggested by 
McClymont & Canfield (1983) is used; AtT > 10^ K, A(T) cal- 
culated by Cook et al. (1989) is adopted. The characteristic val- 
ues taken as basic units for the mass density, temperature, length, 
and magnetic field strength are: po = 3.34 x 10"'" kg m"^, cor- 
responding to the electron number density Ne - 2x 10" cm"-', 
To = W K, Lo = 500 km, and Bo = 25 G, respectively. The 
dimensionless coefficients of the heat conduction and radiative 
losses are in the same forms as given by Ding et al. (2010). 

The dimensionless size of the computational domain is < 
X < 6 and < y < 6, divided into 500 x 500 grid points (in the 
Ding et al. (2010) work a smaller grid of 400x400 was adopted). 
Uniform meshes are adopted in both x— and y- directions. In 
the present study, x is the horizontal axis, and y is vertical, rep- 
resenting the height of the solar atmosphere. At the bottom, all 
quantities are fixed. Other boundaries are treated as open, and all 



quantities are calculated in terms of equivalent extrapolation. A 
multi-step implicit scheme (Hu 1989) is used to solve the MHD 
equations. 

2.2. Initial State 

The initial magnetic field is a potential one taken to be in the 
following form 

\ = 0, 

where if/bo is a free parameter used to control the magnetic field 
strength of the initial background. The initial magnetic field is 
different from that used in Ding et al. (2010) where a linear 
force-free field with a vertical current sheet is adopted. The ini- 
tial state is assumed to be static and isothermal, with a temper- 
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ature T = = 2 x 10"* K representing the atmosphere in the 
chromosphere. The initial mass density is calculated by 



Cose Al at t = 



p = pbo exp(-gy/2). 



(4) 



where pho is a free parameter used to control the mass density at 
the bottom. 

By choosing different values of t(^bo and p^o, we can obtain 
four initial states represented by Al, A2, Bl and B2. The param- 
eters of the four initial states are listed in Table. 1. The initial 
states are divided into two groups, case A (iffho - 2) and case B 
(i^feo = l)j corresponding to the initial backgrounds with strong 
(50 G) and weak (25 G) magnetic field strength, respectively. 
Then, each group is divided into two subgroups according to the 
value of pho- The case of pi,o - 1 (Ne = 2 x 10" cm"-') rep- 
resents the initial state with high mass density, and pbo - Q.I 
(Ne=2x 10'° cm-^) with low mass density. The characteristic 
values of Alfven velocity (va) and the ratio of the gas pressure to 
magnetic pressure (J3) corresponding to each state are also listed 
in Table 1. As the initial temperature of the background is the 
same for each state, the characteristic sound speed for all the 
cases are the same as well, namely 23 km s"'. 

Fig. 1 (left panel) shows the configuration of the initial mag- 
netic field for the case of Al as an example. The initial magnetic 
field is uniform and vertical. Fig. 1 (right panel) shows the dis- 
tribution of the initial mass density versus height for case Al. 
The initial mass density is also uniform along x. As the radiative 
loss is reduced to zero at T = 2 x lO'* K, the initial state is not 
only in hydrostatic equilibrium but also in thermal equilibrium. 

The magnetic configuration of case A2 is the same as case 
Al. The distribution of the initial mass density for case Al and 
A2 are similar, but its value is only one tenth of that for A2 as 
regards the mass density at the same height. For cases in group 
B, the magnetic configurations are the same. The magnetic field 
strength in A is twice of that in B. The distribution of the initial 
mass density of case Bl is the same as case Al, and that of case 
B2 is the same as case A2. 



2.3. New Magnetic Flux Emergence 

The emergence of new magnetic flux is implemented numeri- 
cally in the same way as presented in Ding et al. (2010) (see 
Eq. 10 therein). Here, the flux emergence time is take to be 80 s 
as well. The a parameter which controls the magnitude and ori- 
entation of the emerging flux is taken to be -1.2 for cases in 
group A, and -2.4 for cases in group B. The reason why a val- 
ues are different is to keep the ratio of the newly emerging mag- 
netic flux to the total flux of the initial background the same for 
case A and B, so that reconnection between the new flux and 
pre-existing flux occurs at a similar height in the two cases. The 
localized resistivity introduced to initiate magnetic reconnection 
is in the same form as expressed in Ding et al. (2010). 



2.4. Line profiles 

In order to compare the numerical simulations with observations, 
line profiles are calculated in terms of 



X=4 V= V2 

,0 = J J p-G(T)eT(A)dxdy, 



(7) 



s (see Fig. 2(a)) 

V ^ 100 km 




5 o 
E 



Fig. 3. Re-plot of the state at f = 80 s in Fig. 2(a). The tempera- 
ture is shown in a smaller color scale than in Fig. 2(a). The two 
solid lines (black) denote the .^-position of the slit, and PI and 
P2 denote the two y-position of the slit. 
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Fig. 4. Evolution of the physical quantities, e.g. maximum jet 
velocity (solid line, top row), its line-of-sight component (dotted 
line, top row), and plasma temperature (bottom row) at the loca- 
tion of the peak jet velocity, for case Al . (Left column) measured 
at PI; (Right column) measured at P2. 

where y represents the line-of-sight, G(T) is the emission contri- 
bution function, and 



er(/l) ' 



exp-[^]^ 



(8) 



is the static line profile. In the equation above, Aq is the rest 
wavelength of the resonance line. As - ^Vp is the Doppler shift 
corresponding to the line-of-sight velocity, Vp, and aAq is the 
Doppler width of the line given by 



^0 2kBT 



(9) 



A-=l >'=>'! 



where c is the speed-of-light, kg is the Boltzmann constant, and 
m, is the ion mass. The contribution function, G(T), is obtained 



J. Y. Ding et al.: Magnetic reconnection implications for jet formation 



5 



C III 



V 



Fe IX 



3.5 



2.5 - 



1.5 - 



0.5 




3.5 



2.5 



0.5 



20 s 
40 s 
60 s 
80 s 





150 







150 -150 







150 -150 







150 



3.5 




150 







150-150 150-150 

Doppler Velocity (km s~^) 







150 



Fig. 5. Line profiles for the C iii 977 A line {left column), the O v 629 A line (middle column), and the Fe ix 171 A line (right 
column) at four times for case Al. Calculations are perfomed at two positions; PI (top row) and over the whole region (y - Q - 6, 
bottom row). 



from the ADAS database (Summers 2004), and ionization equi- 
librium is considered. During the calculation of line profiles, a 
1500 km wide slit, located between 1 < x < 4 (see Fig. 3), is 
used. When observed with spectrometers like SUMER or EIS, 
line profiles from different regions of the jet can be obtained, 
e.g. see Figs. 7 & 8 in Madjarska (2011), where the small box 
outlined in Fig. 7 may be moved around to produce line profiles 
in a different part of the jet and/or reconnection site as seen in 
Fig. 8. Therefore, it is important that the simulation derived ob- 
servables from various parts of a jet are studied and compared 
with observed line profiles from the same type of regions. We, 
therefore, divided the simulated regions into two parts: PI which 
covers the region between yi = and y2 = 2, i.e. over the recon- 
nection site and its close surroundings, and P2 between yi = 2 
and y2 = 5, i.e. over the entire jet. The positions of PI and P2 
are also denoted in Fig. 3. Note however, that in some cases, de- 
pending on the jet temperature and density, there may be little or 



no difference in the simulated profiles at the different depths/lo- 
cations, see Al case later Furthermore, if the jet was traveling 
directly towards the observer with no line-of-sight inclination, 
then deriving of the line profiles from different jet depths/loca- 
tions would not be possible. However, jets are seen propagating 
in any directions on the solar disk, including the example out- 
lined above, where such a comparison is possible and therefore, 
needs to be investigated both observationally and theoretically. 



In the present study, we will focus on the jets resulting from 
magnetic reconnection as seen in different parts of the solar at- 
mosphere. Line profiles formed in the lower transition region 
(TR) are calculated for C iii 977 A (8 x 10"* K), the upper tran- 
sition region as seen in O v 629 A (2.5 x 10^ K) and the lower 
corona as viewed in Fe ix 171 A (8 x 10^ K). 
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3. Numerical Results 

In this section, the results of magnetic reconnection triggered by 
newly emerging magnetic flux are presented. Two cases are con- 
sidered, case A and B, where the initial background magnetic 
field strength is different. In each case, two sub-cases with dif- 
ferent initial mass density are investigated. 



Vp - cA JAq. In our line profile calculations the jets propagates 
towards the observer, i.e. the observer is at the top of the y-axis. 
Therefore, blue-shifted emission corresponds to plasma moving 
from the bottom to the top of the y-axis ("+"), while red-shifted 
emission is associated with plasma moving in the opposite di- 
rection ("-"), i-e away from the observer. Whereas in the defi- 
nition of Doppler velocities, blue-shifts are regarded as negative 
values and red-shifts as positive. Therefore, the velocity values 
obtained in our simulations need to be reversed in the calcula- 
tions of line profiles, so as to get the same expression of Doppler 
velocities as in the observations. It is shown that the maximum 
tics iffto = 2 and pho = 1, i.e. B = 50 G and Ne =2x 10" cm"^ blue-shifts reach 60 km s"' at f = 80 s for the C iii 977 A fine. 



3.1. Case A1: cool chromosphere-TR jet 

In case Al we explore an environment with plasma characteris- 



which corresponds to high field strength typical of an active 
region and high density environment. Fig. 2(a) shows the evo- 
lution of the magnetic field, temperature, and velocity for this 
case. As the new magnetic flux with negative polarity {a < 0) 
emerges, the background magnetic field lines will be pushed out 
towards x - Q, resulting in horizontal flows. A current sheet is 
then formed at the right-hand side of the newly emerging mag- 
netic arcade. If magnetic diffusion is introduced into the current 
sheet, magnetic reconnection will occur Part of the magnetic 
field lines expand outward (towards the upper-left in Fig. 2), and 
the plasma at the right-hand side of the flux emergence region 
is pushed into the diffusion region. As reconnection goes on, 
the expansion of magnetic field lines propagates to higher re- 
gions, where plasma is dragged by the magnetic field and moves 
towards the left as seen in the figure. At the same time, mag- 
netic field lines together with plasma from the right-hand bound- 
ary are driven leftwards due to pressure imbalance. This causes 
strong horizontal flows at higher altitude, as seen at f = 60 s and 
80 s in Fig. 2(a). These horizontal flows will collide with the up- 
ward outflows as a result of magnetic reconnection, and a shock 
is formed at the interface. In order to show the shock, the state at 
f = 80 s in Fig. 2(a) is re-plotted in Fig. 3 where the temperature 
is plotted using a smaller colour scale. At the interface of the 
horizontal flows and upward outflows, a bundle of plasma hotter 
than its surrounding is clearly seen. During the magnetic recon- 
nection, the plasma in the diffusion region is heated by Joule 
dissipation. The outflows, as a result of magnetic reconnection 
is bi-directional, upward and downward, along the current con- 
centration. In our model, only part of the down-flow is visible 
because of the low X-point. 

Fig. 4 (top left) shows the maximum jet velocity (solid line) 
and its line-of-sight component (dotted line) measured at PI, as 
a function of time. The corresponding plasma temperature at the 
location of the maximumjet velocity is also shown in Fig. 4 (bot- 
tom left), as a function of time. During the process of magnetic 
reconnection, the jet is accelerated by magnetic tension forces, 
reaching a velocity of 150 km s"' at f = 80 s in the PI re- 
gion, with a line-of-sight velocity of 100 km s"'. The temper- 
ature of the jets with peak velocity is about 4 x 10^ K (20 T,) 
at maximum. These measurements would correspond to a re- 
gion at the foot-points of a chromospheric jet. The same phys- 
ical quantities are also measured at P2, where the jet velocity 
only reaches 45 km s"' at maximum and the jet temperature is 
X 2-3 X 10"* K (1-1.5 r,). These plasma parameters describe a 
short jet reaching only typical chromospheric temperatures prop- 
agating at a speed which is comparable with observed values of 
fibrils. 

Fig. 5 (top row) shows the C iii 977 A line (top left), 
O V 629 A line (top middle), and Fe ix 171 A line profiles (top 
right) measured at PI at four times, t - 20, 40, 60, 80 s. 
The line profiles are plotted as a function of the Doppler ve- 
locity (vp), instead of wavelength (A), in terms of the expression 



and -90 km s ' for both O v 629 A and Fe ix 17 1 A. All the line 
profiles show at least a two-Gaussian structure. 

For C III 977 A, the radiance of the red-shifted component 
increases during the interval f = 20 s to 40 s, and then decreases. 
The red-shifted component centers around velocity 0, which is 
mainly contributed by the plasma in and/or close to the diffusion 
region. From f = 20 s to 40 s, more plasma is heated, but to a 
temperature less than 2 x 10^ K (which is the upper threshold of 
the contribution function for C iii 977 A). As the plasma is heated 
to above this temperature, its contribution to C iii line becomes 
very weak and negligible. The blue-shifted component comes 
from the plasma outside the diffusion region. It is closer to the 
diffusion region at f = 20 s, and has a strong emission. Later, the 
hot region expands, but the high velocity region expands slowly 
and its mass density becomes low, leading to a decrease in the 
blue-shifted component at f = 40 s. As more plasma is heated 
and accelerated, blue-shifted emission starts to increase from 
f = 40 s to 60 s. For O v 629 A, the radiance of the highly blue- 
shifted components increases from f = 20 s to 40 s. This is be- 
cause more plasma is heated to higher transition-region temper- 
atures. As the temperature increases to above 2.5 x 10^ K (i.e. the 
temperature of maximum ionization (T,„ax) for O v 629 A) after 
f = 40 s, the contribution of the hot plasma to the O v 629 A line 



C III 




-150 150 

Doppler Velocity (km s~') 



Fig. 6. Line profiles calculated at P2 for C m 977 A line at two 
times for case Al. The radiance is shown as a function of the 
Doppler shift. 
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10 30 50 70 40 50 60 70 80 

Time (s) Time (s) 

Fig. 7. Evolution of the physical quantities, e.g. maximum jet 
velocity (solid line, top row), its line-of-sight component (dot- 
ted line, top row), and the plasma temperature (bottom row) at 
the location of the peak jet velocity for case A2. (Left column) 
measured at PI, and (right column) measured at P2 . 

becomes weaker Then, the radiance of the highly blue-shifted 
components decreases, which corresponds to the increase of the 
radiance of the blue-shifted components for Fe ix 171 A line. 

Fig. 6 shows the line profiles of C iii 977 A line measured 
at P2, where the blue-shift is 20 km s ' at f = 60 s and reaches 
~30 km s ' at f = 80 s. There is no C ni emission at f = 20 and 
40 s, because the up-flow jets reach P2 region after f = 40 s. 
In addition, when measured at P2, the O v 629 A line and 
Fe IX 171 A have no detectable emission. 

The results in case Al are summarized in Table. 2, where the 
maximum velocity of the jet (Vjet), its line-of-sight component 
(Vy), the plasma temperature (T^^t) at the position of the maxi- 
mum jet velocity, and blue shifts for C iii 977 A (Vc), O v 629 A 
(Vo), and Fe ix 171 A line (Vpe) are calculated at two positions, 
PI and P2, respectively. All the values shown in the table are 
the maximum ones obtained in 80 s. The results for other cases, 
namely A2, Bl and B2, which are discussed in detail in the fol- 
lowing sections, are also listed similarly in Table. 2. 

For the Al case, line profiles are also calculated over the 
whole height, e.g. y = 0-6, as shown in the bottom row in Fig. 5. 
There is very small difference of line profiles at PI and that of the 
whole height. In fact, the line profiles at two different depths are 
the same for O v 629 A line and Fe ix 171 A lines, respectively, 
this is because there is a zero contribution in the O v and Fe ix 
emission above the PI region. 

3.2. Case A2: hot transition region jet 

In the A2 case, the magnetic field background is of the same 
strength as in case Al, i.e. 50 G, but the environment has a 
lower electron density (Ne - 2x 10"^ cm"-'), i.e. typical of a re- 
gion formed higher in the atmosphere than the Al case. Fig. 2(b) 
shows the evolution of the magnetic field, temperature, and the 
velocity. The plasma in the diffusion region is heated to higher 
temperatures, about 1.5 x 10* K (75 Tj) at maximum, compared 
with the case of Al. Because an environment with lower den- 
sity is considered, hotter plasma (i.e. transition-region tempera- 
ture and above) spreads over a larger region (this corresponds to 
about 700 km wide and 3000 km high) at f = 80 s, whereas the 
hotter plasma region is only ~ 100 km wide and ~ 1000 km high 
in the Al fibril case. The decrease in density will increase the ra- 
tio of Joule heating to radiative losses, so that the plasma in the 



diffusion is more strongly heated. Moreover, the hot plasma is 
ejected outward, which will also heat the plasma outside the dif- 
fusion region. As the background atmosphere is at lower density, 
the energy losses of the heated hot plasma outside the diffusion 
region by radiation becomes small, so that the plasma over a 
larger area remains hot. 




10 30 50 70 50 60 70 80 

Time (s) Time (s) 

Fig. 9. Evolution of the physical quantities, e.g. maximum up- 
flow jet velocity (solid line, top row), its line-of-sight component 
(dotted line, top row), and the plasma temperature (bottom row) 
at the location of the peak jet velocity for case Bl. (Left column) 
measured at PI; (right column) measured at P2. 

Fig. 7 shows the maximum jet velocity (solid line, top row), 
its line-of-sight component (dotted line, top row), and the plasma 
temperature at the peak jet velocity (bottom row) as a function 
of time, for the A2 case. At PI (left column), the up-flow jets 
reach 320 km s"' at maximum, line-of-sight component being 
240 km s"'. The temperatures of the up-flowjets with maximum 
velocity reach about 1 .3 x 10'' K (65 T,). It seems that the outflow 
jets with high velocity and high temperature are comparable to 
hot X-ray jets. When observed at P2 (right column), however, 
both the velocities and temperatures are smaller: the jet veloc- 
ity and its line-of-sight component reach about 90 km s"' at 
maximum, its temperature being about 6 x 10^ K (30 T,) at max- 
imum. This means that a jet at transition region temperatures is 
produced. The differences between the jet velocity and its line- 
of-sight component are very small, showing that the outflows 
with maximum velocity are almost vertical. Moreover, there is 
a sharp drop in the temperature of the up-flow jets with peak 

Table 2. Summary of the results for all the cases, namely case 
Al, A2, Bl and B2, where " — " means no signal. See text for 
details. 
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Fig. 8. Line profiles of C iii 977 A line (left column), O v 629 A line (middle column), and Fe ix 171 A line (right column) at four 
times for case A2, where the line radiance is shown as a function of Doppler shift. Each row shows the observations between certain 
heights. Top row: yi = and ^2 = 2 (PI); Bottom row. y\ -2 and ^2 = 5 (P2). 



velocity between t - 15 and 80 s. This is because the outflows 
reach maximum velocity at higher altitudes, where the plasma 
temperature is low. 

Fig. 8 shows the line profiles of the C iii 977 A line (left 
column), the O v 629 A line (middle column), and the Fe ix 171 A 
line (right column) at four times, t - 20, 40, 60, 80 s, as a 
function of Doppler velocity. Line profiles are calculated at two 
positions: PI (top row) and P2 (bottom row). Observed at PI, the 
maximum blue-shifts are of the order of 70 km s ' for all three 
lines. The intensity of C iii 977 A and O v 629 A is a factor of 
2-3 lower than in the Al fibril case, while Fe ix 171 A is larger 
The blue-shifted components of all three lines are stronger when 
observed at P2 than at PI at f = 80 s, whereas the red-shifted 
components is absent at P2. In this case of TR temperatures and 
relatively strong magnetic field, the produced jet is capable at 
reaching high TR temperatures which can peak for a short period 
of time at low coronal temperatures. Depending on the line- 



of-sight, the observer may see a mixture of PI and P2, and in 
instances where the jet is moving directly towards the observer, 
we will see the sum of PI and P2 line profile. 



From an observational point of view, we should expect a 
small-scale brightening to occur at coronal temperatures but a 
jet will only be detected at transition region temperatures. This 
type of phenomena are already seen in multi-instrument obser- 
vations and described by Subramanian et al. (2011, submitted). 
The line profiles composed of multi-component blue-shift and a 
redshift in transition region lines are similar to observable line 
profiles which describe the so-called 'explosive events' (EEs). 
Note, however, that EE spectral line profiles were also associated 
with a surge in the plage area of active region (Madjarska et al. 
2009) as well as were hnked to EUV jets (Madjarska et al. 201 1, 
in preparation). 
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3.3. Case B1: cool chromospheric jet 

Here, the physical environment has a weaker magnetic field 
strength of 25 G, and a high plasma density {Ne - 
2 X lO" cm~^). Fig. 2(c) shows the evolution of the magnetic 
field, temperature, and velocity. The physical process is the same 
as explained in Section 3.1, i.e. horizontal flows and a shock in- 
terface between the inflows and outflows. The plasma in the dif- 
fusion region is heated to 3 x 10^ K (15 Tj). The scale of the hot 
plasma region (i.e. transition-region temperature) is very small, 
50 km in width and 800 km in height, i.e. a sub-arcsec region. 

In Fig. 9 at PI, the velocity of the up-flows {solid line, top 
left) reaches ~100 km s"' at maximum, with a line-of-sight com- 
ponent (dotted line, top left) of ~60 km s"'. The temperature of 
the up-flows with maximum velocity is ~6 x 10^ K (3 T,). At P2 
{right column in Fig. 9), the maximum jet velocity {solid line, 
top right) is only ~40 km s~' at maximum, and the jet tempera- 
ture {bottom right) is 2-3 x 10"* K (1-1.5 T,). Therefore, a cold 
jet that propagate with relatively low speed is produced as com- 
pared to cases Al and A2. 




Fig. 10 shows the line profiles for C m 977 A at four times. At 
PI {left) the blue-shifts reach about 60-70 km s"' at maximum, 
while at P2 {right), the blue-shifts are only about 25-35 km s"'. 
As the hotter plasma region is small in scale, the emission at 
higher altitude is weak or absent. The temperatures is not high 
enough to produce sufficient O v 629 A and Fe ix 171 A line 
emission. A small-scale brightening associated with the diffu- 
sion region should be detected in transition region lines, while a 
small-scale jet could only be seen at low transition region tem- 
peratures at maximum. 



3.4. Case B2: cool chromosphere-TR jet 

The results of case B2 are also analyzed similarly, see Fig. 2(d), 
Fig. 11, Fig. 12, and Fig. 13. In this case, we have a lower 
electron density, A^^, = 2 x 10'° cm"-^, and a field strength as 
in Bl, i.e. 25 G. It shows that the plasma in the diffusion re- 
gion is heated up-to 1 MK (50 Tj) at maximum. Observed at PI 
{left column in Fig. 11), the maximum velocity of the up-flows 
reaches 180 km s"' at maximum, with 130 km s"' in the line- 
of-sight component. The temperature of the up-flow jets with 
maximum velocity reaches 7 x 10^ K {35 Tf). At P2 {right col- 
umn in Fig. 11), both the maximum jet velocity and its line-of- 
sight component are ~60 km s"' at maximum. Moreover, after 
f = 60 s, the temperature of the up-flow jets with maximum 
velocity increases. This is because the outflows get maximum 
velocity at low altitude, where the plasma is heated to higher 
temperature. In this environment, the line profiles of all three 
lines show a Doppler shift of 50 km s ', when calculated at PI 
(Fig. 12). At P2 (Fig. 13), both C iii 977 A fine and O v 629 A 
line show blue-shifts of 50 km s ' as well. The emission in the 
Fe IX 171 A line is absent, however, when observed at P2. We 
should be able to register observationally an event which starts 
with a brightening at coronal temperatures, but jets will only be 
seen at transition region temperatures. 



-100 -50 50 100 -100 -50 50 100 

Doppler Velocity (km s~^) Doppler Velocity (km s~^) 



Fig. 10. Line profiles of C iii 977 A line observed at PI {left) 
and P2 {right) at four times for the case of Bl, where the line 
radiance is shown as a function of Doppler shift. 
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Fig. 11. Evolution of the physical quantities, e.g. maximum up- 
flow jet velocity {solid line, top row), its line-of-sight component 
{dotted line, top row), and the plasma temperature {bottom row) 
at the location of the peak jet velocity for case B2. {Left column) 
measured at PI; {right column) measured at P2. 
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Fig. 13. Line profiles of C m 977 A fine {left) and O v 629 A line 
{right) at four times for case B2 at P2, where the line radiance is 
shown as a function of Doppler shift. 
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Fig. 12. Line profiles of C m 977 A line (left), O v 629 A line (middle), and Fe ix 171 A line (right) at four times for case B2 at PI, 
where the line radiance is shown as a function of Doppler shift. 



4. Conclusions and Discussions 

Numerical studies on magnetic reconnection caused by new 
magnetic flux emergence, occurring at densities typical of the 
chromosphere and transition region, starting in an environment 
with a temperature of 2 x lO'* K (Tj) are presented. The mag- 
netic field strength is set at 25 G and 50 G, which are the typical 
observed values for chomospheric jets like spicules (quiet Sun) 
and fibrils (active region). It is shown that the results are strongly 
dependent on the physical environment: a decrease of one order 
of magnitude of the density leads to a one time increase of the 
maximum jet velocity; an increase of two times the magnitude of 
the magnetic field strength leads to the same velocity enhance- 
ment. Here, the ratio of the amount of the newly emerging mag- 
netic flux to that of the initial background is kept the same in all 
cases considered. If less new flux is emerged, or the flux emer- 
gence time is longer, the maximum velocity and temperature of 
the outflow jets will be weaker Furthermore, both the inclina- 
tion and non-uniformity of the initial background magnetic field 
will have an effect on the velocity and temperature of the out- 
flow. Line profiles of two transition region lines (C iii 977 A and 
O V 629 A), plus the lower coronal line (Fe ix 171 A) are cal- 
culated at two positions: a region which covers the reconnection 
site and its close surroundings, and a region which covers the en- 
tire jet for all four cases. The simulated response in these spectral 
line is an important observable which can directly be compared 
with observed line profiles in various types of jets. 

In the case where we have a lower electron density (A^^, = 
2 X 10'° cm""*), i.e. more typical for the transition region, and a 
stronger magnetic field of 50 G, up-flows with line-of-sight ve- 
locities as high as 200 km s ' and temperatures around 1 MK 
(50 Ti) can be produced around the reconnection site. The jet 
at P2, however, reaches ~90 km s"', but its temperatures is 
only 6 X 10^ K (30 T,) at maximum, i.e. upper transition re- 
gion temperatures. As a result, Fe ix 171 A has weak or no 
detectable emission. The opposite case of high density plasma 



(Ne = 2 X 10" cm ^), typical for the chromosphere and 
weaker magnetic field of 25 G, results in a jet with a velocity 
of -40 km s"' and a temperature -2-3x10^* K (1-1.5 T,), re- 
spectively. Therefore, a weak magnetic field in a high-density 
plasma results in a jet not higher than the typical chromospheric 
values. 

Karpen et al. (1995) studied chromospheric eruptions caused 
by shear-induced magnetic reconnection. This work used a much 
larger box with the shear reaching the X-point after 300 s and 
concluded that the dynamics of chromospheric eruptions were 
dependent on the geometry, shear strength, and local resistiv- 
ity. The present study is an important attempt to test the forma- 
tion of jets originating in the lower solar atmosphere as produced 
only by magnetic reconnection triggered by new flux emergence. 
Here, we do not look into the emergence rate, the emergence 
time or the amount of the newly emerging magnetic flux which 
are important factors that could effect the results of the jets. 
Instead, we study the effect of changing the field strength and 
electron density. The parameters have a significant influence on 
the temperature and velocity of the outflow jets. The results 
demonstrate that magnetic reconnection could be an efficient 
mechanism to drive plasma outflows in the chromosphere and 
transition region. It is shown that if the magnetic field is below 
50 G, the jets can be accelerated to relatively high velocities, but 
their temperatures remain in the typical chromosphere - transi- 
tion region range. This work suggests that intensive studies on 
small-scale jets and the magnetic fields associated with them to- 
gether with a detailed spectroscopic evaluation of the jet plasma 
parameters are needed in order to confirm or reject the driving 
mechanism tested here and permit an evaluation of their contri- 
bution to the mass and energy balance in the solar atmosphere. 

Yokoyama & Shibata (1995, 1996) performed simulations 
of magnetic reconnection based on a 2D magnetic flux emer- 
gence model. A huge amount of new magnetic flux is emerged 
and evolves to the coronal heights where the reconnection with 
the pre-existing background magnetic flux actually occurs. The 
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magnetic field strength and mass density around the diffusion re- 
gion are similar to what we use in case Bl where jets with low 
temperature and low velocity are obtained in our study. In the 
Yokoyama & Shibata study, with the reconnection occurring in 
the corona (note that both heat conduction and radiative effects 
were neglected in their study), hot X-ray jets were produced. 
Furthermore, the emergence of new flux also carried cold and 
dense plasma upward, which will then be accelerated by mag- 
netic tension, generating cool jets. 

Nishizuka et al. (2008) explored MHD simulations of new 
flux emergence, triggering magnetic reconnection by adopt- 
ing similar initial conditions to Yokoyama & Shibata (1995, 
1996), but with a hotter (1 MK) and less dense (10'° cm""*) 
corona. Hot X-ray jets (as high as 10 MK), i.e. hotter than in 
Yokoyama & Shibata (1995, 1996), were obtained. 

Murray et al. (2009) studied oscillatory reconnection by in- 
serting a magnetic flux tube below the solar surface. In their ex- 
periment the background initial magnetic field is ~20 G, and the 
initial density decreases from ~10'^ cm in the lower photo- 
sphere to ~10^ cm"^ in the upper corona. The evolution time of 
the magnetic reconnection is very long (~50 mins) for the recon- 
nection system to get close to an equilibrium state around which 
reconnection reversals occur. Once the system reaches an equi- 
librium, the reconnection ceases. The velocity of the jets shown 
in their paper are slow (~50 km s '). However, as the recon- 
nection occurs in the corona, the temperature of the jets is high, 
reaching a few MK. Heat conduction and radiative losses were 
neglected in their simulations as well. 

In the present study the solar atmosphere is fully ionized. 
In the solar corona this approximation is reasonable. However, 
in the lower solar atmosphere, e.g. photosphere and chromo- 
sphere, the plasma is partially ionised, which means neutrals 
are omnipresent. Due to the collisions with neutrals, in partic- 
ular ion-neutral collisions, a Cowling resistivity is introduced 
into the Ohm's law (Khodachenko et al. 2004; Arber et al. 2007; 
Soler et al. 2009, and references therein). As a result, the Joule 
heating is changed to 2yo„fe = 77;jj+/7,.;5^, where 77, 77,., yy, 7^ are 
the Coulomb and Cowling resistivity, and the components of the 
current density parallel and perpendicular to the magnetic field, 
respectively. As far as a 2D MHD model is concerned, the elec- 
tric currents are perpendicular to the magnetic field, and thus, 
the Joule heating term is simplified to Qjouie = ^cj\ = ^cT- 
In the chromosphere, rjc/rj » 1 (KJiodachenko et al. 2004; 
Leake & Arber 2006), which means that the resistivity is in- 
creased when a partially ionized plasma is considered which will 
result in a lower magnetic Reynolds number Litvinenko & Chae 
(2009) suggests the use of an enhanced value of effective resis- 
tivity, instead of the classical resistivity (//), in studying magnetic 
reconnection in the chromosphere. 

In the calculations of line profiles, ionization equilibrium 
means the contribution function is mainly dependent on the elec- 
tron temperature but weakly on the electron number density. 
This is reasonable for the cases where the velocity of the outflow 
is smafl, e.g. see Peter et al. (2004). Peter et al. (2006) found that 
the assumption of ionization equilibrium is acceptable for cases 
with higher velocity if the temperature gradient is small. In our 
simulations, strong jets are obtained in some environments, e.g. 
case A2. Despite this, the scale of such strong jets with veloc- 
ities greater than 100 km s"' is very small, mainly around the 
diffusion region. Over a much larger scale, the velocities of the 
outflow jet is small, see top left panel (measured at P2) in Fig. 4. 
Furthermore, there is a sharp density increase between the hot 
up-flow jet and the cold background plasma, which is helpful 
for the reduction of ionization and recombination times. Inside 



the hot up-flow region, the temperature gradient is much flat- 
ter and in order to calculate the line emission more precisely, it 
may be necessary to take into account the time-dependent ion- 
ization in further studies, e.g. see work by Sarro et al. (1999); 
Roussev et al. (2001b). 
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